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Abstract: We study the continuum version of Sinai's problem of a random walker 
in a random force field in one dimension. A method of stochastic representations is 
used to represent various probability distributions in this problem (mean probability 
density function and first passage time distributions). This method reproduces already 
known rigorous results and also confirms directly some recent results derived using 
approximation schemes. We demonstrate clearly, in the Sinai scaling regime, that the 
disorder dominates the problem and that the thermal distributions tend to zero-one 
laws. 



August 1998 



1 Current address: IRSAMC, Laboratoire de Physique Quantique, Universite Paul Sabatier, 118 
route de Narbonne, 31062 Toulouse Cedex 



1 Introduction 



Since the introduction of the problem of a random walker in a random uncorrelated 
force field by Sinai ]T[ a great deal of activity has focused on this area. From the 
mathematical point of view, in one dimension, many results may be obtained, for 
example Kesten Q and Golosov || were the first to derive the average value of the 
probability density function. Subsequently, results on the distribution of first passage 
times and consequently the probability distribution of the maximum displacement were 
obtained. The model is one of great interest as it is an analytically accessible test 
ground for the study of non equilibrium phenomena in a non mean field situation 
where activation processes play a clear role. Recently a real space renormalization 
group study has been carried out which leads to results on multitime properties of the 
process [[| , demonstrating explicity the phenomenon of aging in this model. Physically 
the model finds applications in the study of the random field Ising model (which at 
low temperature may be viewed as a system of annihalating random walkers). The 
logarithmic diffusion seen in the Sinai model serves as a paradigm for the dynamics of 
domain walls in droplet models for disordered spin systems ||. A thorough discussion 
of the relevance of the Sinai model in condensed matter physics and a review of the 



early results obtained may be found in |11, 12 



In this paper we analyse the continuum version of the Sinai model, using techniques 



similar to those used in the theory of disordered quantum systems [14| and that have 



already been used in the study of the Sinai model || [nj. We rederive the results already 
known in a very compact fashion which in addition sheds more light on the nature of 
these solutions. We derive explicit representations for the distributions (and not just 
the disordered averaged values). We shall see for example that the first passage time 
distributions become zero-one laws, i.e. the thermal fluctuations become insignificant 
and everything is dominated by the disorder. 

The version of the model we shall study is that of a diffusion X t in a random force 
field. The equation of motion is the Langevin equation 

?£ = T H + 4t(x t ), (l) 

where rj is a thermal white noise in time, i.e. it is Gaussian of zero mean with correlation 
function (rjtrjt') = 2TS(t — t'). The force <ft' is derived from a potential <fr(x) which is a 
Brownian motion of zero mean such that (fi'(x)(j)'(y) = a 2 5(x — y). The result that the 
diffusion is logarithmic may be obtained through the following well known argument: 
The time to move from the origin to a point x is given by the Arrhenius law, i.e. 

t ~ exp((0(x) - 0(O)/T) ~ exp(av^/T), (2) 

by Brownian scaling, and hence X{t) ~ log 2 (t). Exact results derived show that this, 
rather crude, argument in fact reflects very accurately the behaviour of the diffusion 
in the appropriate scaling regime. In this paper the logarithmic nature of the diffusion 
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is derived from first principles without recourse to the Arrhenius law or any scaling 
arguments. 



2 Mathematical Formalism 

We start with the model ([!]) of a particle moving in a quenched random white noise 
force field 4>'(x) such that 4>'(x)(p'(y) = 25(x — y), here the overline indicates the average 
over the disorder (i.e. without loss of generality we have taken a = \/2). The particle is 
also subject to a random white noise in time such that the probability density function 
for a particle, started at the origin at time zero, to be at x at time t is given by the 
standard forward Fokker-Planck equation: 

d_ ( dp + d4 \ = dp 
dx \dx dx J dt 

In what follows we shall find it useful to work with the Laplace transform of the density 
with respect to time which is given by 

where E is the Laplace transform variable congugate to t. The adjoint p(x) of the 
probability density p(x) (the probability to be at given that the particle starts at x) 
is related to p(x) via detailed balance 

p(x) = p(x) exp (— (j){x) + 0(0)) (5) 
and can be shown to satisfy the backward equation 

d 2 p d(f>dp , , ( 

— -—— = Ep- 5{x). (6) 
dx z dx dx 

We now proceed by making essentially a Hopf-Cole transformation of the problem 

p(x) = Aexp(—J z^ds^j for x > 

p{x) = Aexpf—J z~ds^j for x < 0, (7) 
where it is easy to see that 

^ = ~E + z ±2 {s) + ^^zHs) s > 0. (8) 
ds ds 

The jump at x = fixes the value of A to be 
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Zq "f Zq 

Hence we have found an explicit representation for p (and hence p) in terms of the two 
stochastic processes z ± (s): 

p(x) = -zr — exp (— [ ztds] for x > (10) 

Zq + Zq \ JO J 

pix) = -jt~ exp f — / z~ds] for x < (11) 

W ZjJ" + Zq \ JO 8 J ~ ^ ' 

In general one requires boundary conditions for z , here we imagine that the diffusion 
is restricted to a finite segment of the real line [-L, L] with reflecting boundaries at 
x = ±L, hence 

^|±l = 0^ zt = zl = 0. (12) 

On using the equation of detailed balance ([5]) and the integrated form of @, one may 
show 



p(x) = + u — exp I — / — ds J for x > 




o zl 



p(x) = + u — ■ exp I — / — ds ] for x < 0. (13) 



z -x\ z + 2 o 

A useful property of p is that we may write 



P(*) = -4 , + f j-expf- (14) 




(2q + 2 ) fix 

using this relation it is straightforward to explicity verify that 

p{x)dx = —, (15) 

i.e. the conservation of probability. The observant reader will notice that the stochastic 
processes z^ have the pathology that they have no normalisable equilibrium measure 
and in addition that they may become infinite (finite time explosion). Let us remark 
here that, for probabilistic reasons, z must remain positive as p(x) must decrease 
monotonically on increasing the distance from the origin. In fact the natural process 
to work with is wf = zf_ s , which obeys the time reversed (with respect to z) stochastic 
differential equation 

dw 2 dxf> 

— = E — w+—w. (16) 
ds ds 
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Using the Stratonovich prescription the generator for the diffusion w is 



G w = ^-w^-w-(w 2 -E)^-. (17) 
ow ow ow 



The equilibrium distribution for w therefore satisfies 

QlP (w) = 0, 

and is hence given by 



PoH = 1 , -exp(-w - -), (19) 
K (2E2)w w 



where Kq (x) is a Bessel function of the third kind ]T5| . It may be convenient , depending 



on the situation, to work with the alternative variable 

w = £5exp(7), (20) 
the fact that we are using the Stratonovich prescription means that we may apply the 



chain rule to the stochastic differential equation (16) and hence find: 

^ = -2£§sinh( T )-^ (21) 
as as 

Slightly abusing notation, the equilibrium distribution for 7 is given by: 

^o(7) = — ^tt exp(-2£* cosh( 7 ))). (22) 
The boundary conditions induced on w are: 

^|±l = -> zt = wt = (23) 
For example, for x > 0, p is given in terms of w as 



p(x) = + _ exp - / wfds . (24) 

In the limit L — > 00 the boundary conditions at ±L should become unimportant. The 
following subtle argument demonstrates this fact. On sending L to infinity (whilst 
keeping x finite) the distribution of wl is given by the equilibrium measure (|19|) (i.e. 
Wl is in equilibrium because L corresponds to a large time by which the process w is 
in equilibrium). In addition, as L — > 00 the distribution of wl- x is also given by the 
equilibrium distribution. Hence treating L — x as the new origin for the process w + , L 
as a new origin for w~, and using the independence of w + and w~ we may write the 
statistically identical expression: 

p(x) = — — — - exp (— [ wjds) , (25) 



wj + w V Jo 

± 



where it is to be understood that w are in equilibrium. 
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3 Probability Density Functions 



With the formalism elaborated in the previous section it is now immediate to calculate 
the probability density at the origin. This result has been calculated using the same 
technique as here in |3|, [Kj and we include it here for the sake of completeness as its 
derivation is now immediate. The same result may also be found in |13|| where it is 
calculated using the replica method from disordered systems. 

- P {w)P {w') 



p(0) = / uv ' uv ' dwdw', (26) 
v ' Jo w + w' v ' 

and we start by expressing the denominator in the integral (p6| ) as the integral over an 
exponential, i.e. 



P(0) = -r— ^— t- r exp (-t(x + y )-E 1 Hx + -)-EHy + -)) ^-dt. (27) 
E2K (2E2) 2 Jo \ x y J xy 

The integrals over x and y may now be performed yielding (in terms of Bessel functions 
of the third kind |]J) 

1 ' --2/ 



P(0) = -T— — T7T / h Kt( X )xdx. (2J 



We now use the differential identity 



U^ { Kl-Kl))= X Kl (29) 



thus yielding the result: 



AKA2E 



h\2 



p(0) = 1V i -4. (30) 

K Q {2E2) 2 



Using the asymptotic expansion for Bessel functions of the third kind [15] we find that 
for small E (corresponding to large time) 



" 0) ~ B*m< (31) 

and hence using the Tauberian theorems for regularly varying functions we obtain the 
result 

1 



p(0,0,*)~— 5777 . (32) 
log (t) 



To examine the full (x, t) dependence we note that 
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^o(7o + )^o(7o" 



ft?) = / T^i ; + z, 1 } ^ G ™> 7 ' W^7o~, (33) 

•/ E'2 exp(7+) + E* exp(7 ) 

where we have used the Feynman-Kac Formula and G is the Greens function obeying 

dG d 2 G i , , . dG i . . _ , . 

2^2 sinh(7o)^ £2 e xp(7 )G, (34) 



<9x (97q c?7o 
subject to the initial condition 

G( 7o ,7,0) = <S(7-7o)- (35) 

In general the analysis of this Greens function is difficult and few analytic results 
are known. However if we restrict ourselves to the regime of Sinai scaling i.e we 
write x = X\og 2 (E) and 7 = —u\og(E), where u and X are of order one, we find a 
considerable simplification occurs which permits a rather thorough asymptotic analysis 
of this regime; this is the key remark that enables us to extend the method of [[J to 
obtain new results. In these new variables, the equation ( 34) becomes: 

art FfiG r)G 

— = ^± + - E*+«) M£)|^ - \og 2 (E)E*- u °G. (36) 

In the region where w > | we find the equation is dominated by the potential and 
drift terms i.e. 

H = \og{E)E-^^- - \og 2 (E)E^°G, (37) 



which has the general solution 

^^w" 1 ^ 1 ' (38) 

and hence G(uq,X) ~ for uq > ~. In the region uq < — |, the equation is dominated 
by the drift term, i.e. 

the general solution to this equation being G(uq,X) = f(X H — - — ) hence 

G(uo, X) ~ f(X), i.e. it becomes independent of u . In the interval [— ~, ~] the 
equation is dominated by the diffusive term, i.e. 

dx = M- m 

If we apply detailed balance in (|33"D i.e. use 
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n(7o + )G(7o + , 7 + ; x) = P (7 + )G( 7 + , 7o + ; x) (41) 

then we find 



p(x) = / — uw 7 u \ ,u/ -F{ 1 + ,x)d 1 + d% (42) 

•/ exp(7+) + E2 exp(7 ) 

where _ 

F( 1 + ,x)= G( 1 + , 1 +,X)d 1 +. (43) 



By the Feynman-Kac formula F(j,x) obeys 

OF d 2 F 1 , / \ 9F 1 

— = — -2E* smh( 7 )— - £5 exp( 7 )F (44) 

with induced initial condition F( 7 , 0) = 1. Returning to the variable w we finally 
obtain: 

°° exp(-w -w'-R-R) w 



p(x) = — / — - 7 ^ — ^F(log(— ),x)dwdw'. (45) 

K (2E^) 2 Jo ww>(w + w') v &v £^ 

We now exploit the fact that E is small 



1 r°° exp(—Ez — Ez' — - — h) 



p{x) = — / — ^ z -^F(\og(zE2),x)dzdz' 

EK (2EhyJo zz>(z + z>) v 6V ; ' ; 

1 r 00 exp' 1 ' 



/ — - ^—F(— 00, x)dzdz' 

Jo zz'(z + z') 

exp( — u — v! — a(u + u'))dudu' da 



EK Q (2E2yJo zz'(z + z 
F(—oo, x) 
EK (2E^) 2 
F(— 00, x) 
EKq{2EI\ 
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(46) 



From our previous analysis we find that for u G [— |, ~] 

<9F <9 2 F 



«9X 9m 2 ' ^ 
For Mo < — |, -F is constant and continuity at u = — | implies that ^| u= _i = 0; hence 
F(— 00, x) = F(—~,x). Continuity at u = ~ yields = 0. The equation (|47|) for 

F in [ — 95 5] ma y ^ e easu y solved in terms of Fourier modes to yield 

x 4^ 1 /2n + 1 1 \ , (2n+l) 2 „ x 

F(u,X) = - £ — an ((— )( 5 - »)») exp(-L_X^). (48) 
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Note that if one had a Brownian motion in [— ~, ~] started at a position u with an 
absorbing boundary at u — \ and a reflecting boundary at u = — |, then F(w,X) is 
simply the probability that the particle has not been absorbed before time X, indeed 
the function F appears directly in the solution of Kesten Q with this very interpretation 

In the Laplace transformed variables one finds that 

— 4 ~ (-1)™ . (2n + l) 2 2 x , 

7r£/ log (E) ^ 2n + 1 4 log 

Putting in the value of F in equation ( pf) and asymptotically inverting the Laplace 
transform using the Tauberian theorems we obtain the result of Kesten [|J: 

_ 4 ~ (-l) n , (2n + l) 2 „ x , 

Trlog (t) n=0 2n + 1 4 log (t) 



4 First Passage Time Distributions 

In this section we investigate the distribution of first passage times for the Sinai dif- 
fusion, these results of course give information of the distribution of the maximum of 
the process. 

Let us define P(x, t) to be the probability starting at x that the diffusion reaches the 
origin before time t. It is convenient to work with P evaluated at a random exponential 
time Te of rate E, i.e. 

P{x,E)= P(x,t)Eexp(-Et)dt= lexp(-ET(x))), (51) 
Jo 

where the angled brackets indicate the thermal average (i.e. the average over the white 
noise rj) and T(x) = inf{t : X t = | X = x}, i.e. the first passage time to starting 
from x. A recurrence equation may now be found for P(x, E) as follows, starting the 
process at x evolve the particle for a time dt, then 



P(x,E) = ((1 - Edt)P(x + dX t ,E)). (52) 

Here the first term is simply the probability that the time Te does not occur during the 
first step and the second term is P evaluated at the new position x + dX t . Expanding 
P(x + dXt, E) to order dt and averaging over dX t , one obtains the first passage time 
equation 

dx 2 dx dx ^ 
Using the boundary condition P(0, E) = 1 and making the Hopf-Cole transformation 
of the previous section we find: 



8 



P{x, E) = exp(- r z s ds). (54) 
Jo 

where z s obeys the stochastic differential equation @. Re-expressing the above in 
terms of the process w we find 

P(x,E) = exp(— / w s ds). (55) 

J L—x 

Consider x = L, i.e. we consider the time taken to cross the sample with reffecting 
initial conditions at the starting point (hence zl = Wq = 0). If we consider average 
value of the q th power of P(L, E) i.e. P q (L, E) for fixed L and E, using the Feynman- 
Kac formula one finds that P q (L, E) = F q (L, E) where 

8F d 2 F i BF i 

-q£ = -gj - ™* siuh( 7o )^ - q E* exp( 7o )F 9 . (56) 

It is now important to note that for q > 0, in the Sinai scaling regime, the prefactor 
q does not alter the solution for F q , hence one has that P q (L,E) = P(L,E), that is 
to say that P tends to either zero or one independently of the thermal noise i.e. the 
distribution of P(L, E) for fixed L and E, V tends to the form V(P) = aS(P — 1) + 
(1 — a)S(P). Using the calculation of the previous section for F we find therefore, in 
the regime of large times and in the Sinai scaling regime where L/log 2 (£) -0(1), 
that 

4 ~ (-If , (2n + l) 2 , L N 

a = -Y- J —exp(-± —^n 2 — _-). 57 

7r^ 2n + l U 4 log 2 ^)^ 1 J 

One may also consider the case where L — > oo, in this limit we see that 

P(x, E) = exp(- [ X w s ds). (58) 
Jo 

but where the boundary condition on wq is that it is in equilibrium with the measure 
(P^|), here we therefore find 

/OO rp 
Poho)F( lo ,—^—). (59) 
-oo log (E) 

In the Sinai scaling regime xj log 2 (-E) ~ 0(1) we may therefore write: 

I /-oo , , 

P(x,E) = — — - \ exp(-E5-«o - E-^)F(u )\log(E)\du , (60) 
K (2E2) j-oo 

which becomes, in the limit of small E, 

— f\ , , , 8 ~ 1 / (2n+l) 2 A \ 
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Here it is interesting to note that one may calculate explicitly log(P(x,E)) which is 
given as 

\og(P(x,E)) = -x / P (w)wdw , (62) 

J 

which in the limit of small E (or large time), but not dependent on being in the Sinai 
scaling regime, yields 

hence, whereas the scaling of the integer moment of P(x, E) was as y og i^ ; the scaling 
of the average of the logarithm ggesting that in a typical sample the 

scaling should be as log ^ ■ If we consider the case where x is infinitesimally close to 
the origin we see that the probability to return to the origin (simply developing the 
exponential in equation fl58|) to first order in x) is 



x 



P ^ =1 -Mir (64) 

and hence following the definition of || the probablity of not arriving at the origin 
behaves as p pr (L) ~ L~ 9 where L is the length scale L = log 2 (E) and 9 is therefore 
the persistence exponent for the diffusion and is given by 9 = |, thus confirming the 
results of g | . 



5 The Mean First Passage Time 

In the Sinai model we recall that there are two kinds of average, a thermal average over 
the realisations of the temporal white noise and also the average over the realisations 
of the disorder. There has been much interest in the distribution of the mean first 
passage time (MFPT) . Various authors have remarked upon the multifracticality of its 
distribution, in this section we present an extremely simple argument to demonstrate 
this fact. Consider the case where there is a reflecting boundary at L, using the 
formalism of the previous sections or using the formula for MFPTs in we have 



f- 2 (T L (x))-^(T L (x)) = -l, (65) 
dx z ax ax 

where the subscript L indicates the reflecting boundary conditions at L and Tl(x) is 
the first time the particle arrives at given that it starts at x (the angled brackets 
indicating the thermal average). The solution to this equation for a particle starting 
at L is 

(T L (L)) = [ dx exp(0(x)) f dyexp(-(j)(y)). (66) 

JO Jx 

The moments averaged over the disorder are 
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(T L (L))i = f J] dx ndVn exp f E( min (^> x i) ~ m HVh Vi) ~ 2 minfc, Vj )) ) 

J{x n <y„} 1 \ ij J 



I n 

J{x n <y n } 1 



dx n dy n exp I + If* ~~ ~ 2 N ] • ( 67 ) 
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Hence we see that the moment (T^x)) 9 is the partition function for a one dimensional 
anti-Coulomb gas (in the sense that like charges attract and opposite charges repel), 
where each particle has its antiparticle which is conditioned to be on its right. Rescaling 
the spatial variables we obtain 



r q ( L 

(T L (L))i = L 2n / Yl dx n dy n exp -- ^fla;* - Xj\ + \y t - yj\ -2\x t -yj\) 

J{x n <y n } i \ * 



L , 
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(68) 

where all the integrations are now in [0,1]. When L becomes large this corresponds 
to the limit of low temperature and the partition function is dominated by the lowest 
energy state which is clearly when x n = and y n — 1 for all n, i.e. all the like charges 
condense together to form two macrocharges which repel and hence lie on the two 
boundary points of the unit interval. This energy is simply q 2 L. Hence 



log{(T L (L))*) ~q 2 L, (69) 

clearly demonstrating the multifracticality. The distribution is infact approximately 
log- normal. 



6 Probabilistic Interpretation of Results 

Here we shall interpret the results derived in this paper in a probabilistic fashion which 
gives us some intuition for the structure of the Sinai problem. In terms of stochastic 
processes we shall consider the process u x = —^ x j \og(E); if we define the rescaled time 
variable x = X\og 2 (E) then the stochastic equation obeyed by u is 



Using Brownian scaling we may write 0(Xlog 2 (£)) = <f)(X)\log(E)\ and hence the 
induced stochastic differential equation on u as a function of the rescaled time X is 

^ = - log(E)(E^ - E^) + A 0(X ) (71) 
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In the region outside the interval [—5, 5] the process is dominated by a drift which 
returns the particle to the interval. In the interval [— |, |] the process is simply a 
Brownian motion. If we now consider the functional critical to our study of this problem 

H{X) = exp (- jT w a dsj = exp ^- \og 2 (E) J* E^ Us ds \ , (72) 

we see that as E — > 

H(X) — > 1 if max(u s ) < ^ 

— > if max(-u s ) > -, (73) 

thus clarifying at an intuitive level the zero-one type results obtained in the previous 
sections. 



7 Conclusions and Perspectives 



We have shown that the method presented in |3|, [K| (and similar to those used in the 



Quantum Mechanics of one dimensional disordered systems |f[4}| ) may be extended to 
produce explicit results on spatio-temporal probability distribution functions in the 
continuum version of the Sinai model; notably the result of Kesten for the averaged 
probability density function and also for various first passage time probabilities. The 
method demonstrates the preponderence of zero-one type laws where the disorder dom- 
inates. The continuum formulation of the problem also allows a very simple demon- 
stration of the multi-fractal nature of the mean first passage time. 
The great advantage of the technique is that it allows one to represent transition 
probabilty densities as functionals of a well defined stochastic process and therefore 
quantities depending on several times may in principle be analysed, thus permitting 
the study of aging in these systems. Work on these extensions is currently under 
progress JlS| . 
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